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Abstract 

Previous theoretical work on the boundary layer receptivity problem has utilized 
large Reynolds number asymptotic theories, thus being limited to a narrow part of 
the frequency - Reynolds number domain. We present an alternative approach for 
the prediction of localized instability generation which has a general applicability, 
and also accounts for finite Reynolds number effects. This approach is illustrated 
for the case of Tollmien-Schlichting wave generation in a Blasius boundary layer due 
to the interaction of a free-stream acoustic wave with a region of short-scale varia- 
tion in the surface boundary condition. The specific types of wall inhomogeneities 
examined are: regions of short scale variations in wall suction, wall admittance and 
wall geometry (roughness). Extensive comparison is made between the results of the 
finite Reynolds number approach and previous asymptotic predictions. This com- 
parison also suggests an alternative way of utilizing the latter at Reynolds numbers 
of interest in practice. 


1. INTRODUCTION 

Laminar- turbulent transition is a result of instability of the laminar state. At high 
Reynolds numbers, suitable perturbations to this state are amplified, eventually 
leading to the stable (in-the-large) turbulent state. It is important to note that the 
introduction of these “suitable” perturbations (i.e., instabilities) into the boundary 
layer flow is a necessary prerequisite for transition to occur. The process by which 
the boundary layer internalizes the external disturbances in the form of instability 
waves is known as the boundary layer receptivity. Once generated, these instabil- 
ity waves undergo linear amplification and nonlinear interactions before the flow 
becomes fully turbulent. 

Therefore, receptivity signifies the genesis of boundary layer transition. Goldstein^ 1 ,2 1 , 
Rubanl 1 * 3 ! and Goldstein et afM utilized high Reynolds number asymptotic methods 
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to elucidate the basic mechanisms responsible for making the boundary layer “recep- 
tive” to external disturbances. This pioneering work provided much of the impetus 
for further work in receptivity during the past few years. In the remainder of this 
section, we first briefly describe the general features of Goldstein’s receptivity the- 
ory. Following this, we summarize the receptivity research that followed Goldstein’s 
work, and present the problem addressed in this paper. 

The classical linear stability theory is based upon the disparity between the stream- 
wise length-scales of the instability motion and the mean boundary layer. The 
instability waves correspond to the short-scale, nearly periodic eigensolutions of the 
slowly developing boundary layer which have a small growth rate for parameter val- 
ues between the neutral boundaries. The obvious implication of the locally periodic 
nature of the instability wave is that the wave is decoupled from any other spatially 
periodic motion with a different wavelength. 

It is obvious that in order to generate an instability wave of a particular frequency, 
one would require a forcing which has not only the same frequency, but also a 
spatial scale that matches the instability wavelength. However, free-stream distur- 
bances are governed by the inviscid dynamics outside the boundary layer, while the 
instability waves represent free oscillations (or eigenmodes) of the flow within the 
boundary layer. Hence, the wavelengths of these two types of disturbances are, in 
general, quite different. In fact, in low-speed flows, the free-stream disturbances 
(which can either be acoustic waves or converted vorticity disturbances, i.e., free- 
stream turbulence) have wavelengths which are much longer than the instability 
wavelengths. The free-stream unsteadiness can then be “tuned” to the instabil- 
ity length-scale only through an interaction with the spatial spectrum of the mean 
boundary layer. 

Since a slowly developing (or “quasi-parallel”) mean boundary layer does not have 
the short length-scales required for the tuning process, receptivity will usually oc- 
cur only in regions of non-parallel mean flow. Goldstein showed that there are 
two classes of regions where the mean flow becomes non-parallel. The first class of 
non-parallel mean flow regions corresponds to the region close to the leading edge, 
where the boundary layer thickness is changing rapidly. The second category, which 
is much more diverse, involves regions farther downstream of the leading edge. Here 
the mean-flow becomes non-parallel due to either a short-scale variation in the sur- 
face boundary condition or an adverse pressure gradient provoking a separation of 
the mean boundary layer. The short-scale variations in surface boundary condition 
correspond to changes in surface geometry (wall roughness), wall suction/blowing 
velocity or wall temperature. The latter two cases are especially relevant to laminar 
flow control (LFC) where suction and heating (water) or cooling (air) are utilized 
for stabilizing the boundary layer. The variations in boundary conditions can be 
either local (isolated roughness elements, well-separated suction/heating strips) or 
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distributed over a large number of instability wavelengths (distributed roughness, 
surface wavincss, closely spaced suction/heating strips); however, the role of the 
short-scale variation is the same in both cases, viz., to provide the tuning required 
for producing the appropriate forcing. 

In comparing the two classes of receptivity regions, it is important to note that 
instability waves generated near the leading edge undergo an exponential decay 
upstream of the lower branch neutral stability point, even a weak receptivity mech- 
anism close to the neutral point is likely to be more important than a much stronger 
receptivity mechanism near the leading edgel 5 !. j n practice, the exponential decay 
between the leading edge and the neutral stability point is decreased somewhat due 
to the presence of an adverse pressure gradient region close to the leading edge. 
However, in this paper we concentrate on the second class of receptivity regions. 

Goldstein^ 2 ! and Ruban^ 3 ! presented analytical solutions for the specific case of recep- 
tivity due to the interaction of a free-stream acoustic wave with a small but sudden 
variation in the surface geometry. Their analysis was limited to two-dimensional 
flows in terms of both the mean boundary layer as well as the instability wave 
generated due to the interaction. However, the asymptotic framework of Gold- 
stein and Ruban was later extended by other investigators to a variety of other 
problems in the second class of non-parallel mean flow regions. Bodonyi et 
considered the case where the mean flow perturbation due to the wall hump is not 
small enough to permit linearization with respect to the oncoming boundary layer. 
Choudhari^ analyzed receptivity mechanisms arising in LFC applications which 
employ suction through a porous surface as a means of stabilizing the boundary 
layer. He showed that in addition to the receptivity caused by mean flow gradients 
due to variations in the wall suction distribution, the short-scale variation in the 
admittance of the porous surface directly scatters energy from the acoustic wave to 
the instability wave. The latter mechanism is operational even in the absence of a 
mean flow adjustment. Choudhari and Keischen^ studied the three-dimensional 
interaction of a free-stream acoustic wave with a wall inhomogeneity of above types 
which has short-scale variations in both the streamwise and spanwise directions. 
KersclienW and Choudhari and KerscherJ 9 ! investigated the generation of instabil- 
ity waves near a wall hump due to a convected vorticity disturbance in the free 
stream. Finally, localized generation of T-S waves in two and three-dimensional 
supersonic boundary layer's was analyzed by the present authors^ 10 ! in the context 
of two and three-dimensional supersonic boundary layers. 

All previous investigations discussed above utilized the asymptotic (triple deck) 
framework. Although useful in their own right, these predictions were based on 
a single-term asymptotic expansion. However, the critical receptivity region cor- 
responds to locations upstream of or close to the lower branch of the neutral sta- 
bility curve, i.e., Reynolds numbers smaller than those in the main part of the 
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unstable regime, where the quantitative accuracy of the asymptotic predictions is 
questionable. In this paper, we present an alternative approach for the prediction 
of receptivity in localized regions which would yield relatively accurate results in 
the Reynolds-number-range of interest in practice. For simplicity, attention is re- 
stricted to the case of incompressible flows over two-dimensional airfoils. However, 
the present approach also provides a viable framework for the prediction of localized 
receptivity in such complex flows as compressible and three-dimensional boundary 
layers. In section 2, we formulate the problem for the local interaction of a free- 
stream acoustic wave with an arbitrary wall irihomogeneity. It is shown that the 
local receptivity problem reduces to an inhomogeneous Orr-Sommerfeld (henceforth 
0-S) problem in the Fourier transform space, and the amplitude of the generated 
instability wave can be determined as the residue of the pole corresponding to 
that particular eigenvalue of the 0-S problem which represents the wavenumber 
of the Tollmien-Schlichting (T-S) instability wave. In section 3, we present results 
which show the variation of the receptivity coefficient in the frequency parameter - 
Reynolds number space. Detailed comparison of these results with the asymptotic 
predictions is also presented in Section 3. Although the discussion in this paper 
pertains only to localized inhomogeneities, in section 4, we also point out the pos- 
sible application of this analysis to distributed regions of short-scale variations. In 
addition to accounting for the finite Reynolds number effects, the present approach 
has other advantages over the asymptotic theory, and these are discussed in Section 
5. 

2. ANALYSIS 

This section provides a summary of the finite Reynolds number approach. In order 
to assess the finite Reynolds number effects on localized receptivity mechanisms 
in the simplest possible setting, we focus attention on a two-dimensional, incom- 
pressible flow past a semi-infinite, flat-plate airfoil. However, the present approach 
can easily be applied to other types of geometries as well. The free-stream speed, 
density and kinematic viscosity corresponding to the oncoming flow are denoted by 
U p* and v * , respectively. The unsteady free-stream disturbance is assumed to 
be a plane, harmonic acoustic wave propagating parallel to the plate in the down- 
stream direction. Since the acoustic wavelength is infinite in the zero Mach number 
limit, the free-stream flow consists of a uniform flow with mean velocity plus 
a harmonic perturbation of amplitude u* iC and frequency u>*. The uniform velocity 
fluctuation u* ac is also accompanied by a uniform, time harmonic pressure gradi- 
ent with an amplitude We assume the nondimensional amplitude of the 

free-stream fluctuation to be sufficiently small ( ej s = u * a JV ^ << 1) so that the 
unsteady motion can be treated as a small perturbation of the local mean flow. 
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The receptivity is assumed to occur due to a localized wall inhomogeneity involving 
a variation in one or more of the surface boundary conditions. Some remarks 
concerning the application of results for the localized case to an extended region 
of receptivity are presented in Section 4. We assume that the distribution of the 
normal component of the steady or unsteady velocity at the wall, or the wall height 
distribution has a short-scale variation with respect to the streamwise coordinate 
in a local region which is a distance l* from the leading edge (Figs, la-lc and Fig. 
2). While the nonzero mean normal velocity represents the wall suction/blowing 
distribution, the short-scale distribution of the unsteady normal velocity at the 
wall has been introduced in order to model the effects of rapidly varying surface 
admittance in an incompressible flow. 

The admittance of a surface is defined as the ratio, at the surface, of the unsteady 
normal velocity to the fluctuating component of pressure. Since the absolute value 
of the pressure becomes irrelevant in the incompressible limit, we specify the dis- 
tribution of the unsteady normal velocity at the wall instead of directly imposing 
a short-scale distribution of the wall admittance. Even in the compressible case 
analyzed by Choudhari^ using asymptotic theory, the mathematical statement of 
the inhomogeneous boundary condition reduces to the specification of the unsteady 
normal velocity at the wall. Therefore, the present treatment of the wall admit- 
tance variation poses little difficulty in the comparison of present results with the 
asymptotic predictions. It is worth noting that in addition to modelling the wall ad- 
mittance variation, the unsteady mass flux at the wall may also represent unsteady 
disturbances within the sxiction system, especially the designs with large suction 
slots. 

The maximum amplitudes of the local variations in different types of wall bound- 
ary conditions (as described above) are denoted by the nondimensional parameters 
(w , where the index j varies from one to three depending on the type of wall in- 
homogeneity. Specifically, the parameter <44 corresponds to the normalized wall 
suction/blowing velocity, while e„4 denotes e*/u* c , the unsteady vertical 

velocity at the porous surface normalized by the free-stream disturbance velocity. 
The parameter 44 can be related to the amplitude of the wall- admittance varia- 
tion, where ft*, denotes the wall admittance as defined in the preceding 

paragraph. Finally, 44 denotes the roughness height, H*,/6 *, where 6* is the 
displacement thickness of the Blasius boundary layer at the location of the wall 
inhomogeneity (i.c., the mean flow in the absence of any localized variations). All 
three parameters 44 , arc assumed to be small enough to allow linear analyses of the 
different mechanisms and superposition of the results. While the assumed orders 
of magnitudes for the wall-suction and wall-admittance are consistent with typical 
parameter values in LFC systems, the roughness heights can vary over a wide range 
in practice. The assumption of small roughness height allows us to compare results 
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with the asymptotic predictions of Goldsteim-b 

Since the roughness clement and the variation in wall suction are assumed to pro- 
duce only a small perturbation to the Blasius boundary layer the stability proper- 
ties of the perturbed mean flow (i.e., the mean flow in presence of the perturbation 
due to the wall suction or wall geometry variation) are identical to those of the 
Blasius boundary layer, to the leading order. We assume the Reynolds number 
Jl 6 , = U^6*/u*, based upon the local displacement thickness of the unperturbed 
mean boundary layer, to be large enough such that the stability of the unperturbed 
mean flow at the location of the wall inhoinogeneity r is governed by the classical, 
quasi-parallel stability theory. Strictly speaking, the large Reynolds number as- 
sumption is inherent within the boundary layer approximation for the mean flow. 
In practice, however, the quasi-parallel stability theory has been found to be rea- 
sonable only for Reynolds numbers much larger than those at which the mean flow 
is accurately described by the boundary layer theory. For the Blasius boundary 
layer, this constraint (approximately) corresponds to Reynolds numbers larger than 
the minimum critical Reynolds number of 520. 

We further assume that the length-scale of the local variations, L*, is of the same 
order as the local T-S wavelength. This provides the necessary ingredient for the 
wavelength reduction from the free-stream disturbance to the instability wave. In 
this paper we restrict attention to a localized region of receptivity'. However, remarks 
concerning the extension of these results to distributed regions of receptivity will 
also be presented in Section 4. 

Since the instability' wavelength is much shorter than the length-scale of the mean 
boundary' layer (i.e., f( = L* j (* << 1), we can utilize the method of matched 
asymptotic expansions to simplify the problem. Thus, we match a local approxima- 
tion valid in the region of receptivity to an outer (or “global”) solution valid farther 
downstream. The governing equations in the outer region are homogeneous; hence 
the far downstream solution for the short-scale, unsteady motion produced by the 
interaction of the free-stream disturbance with the wall inhomogeneity corresponds 
to a superposition of the discrete and continuous spectrum ty'pe eigenmodes in a 
slowly developing boundary layer. The amplitudes of these eigensolutions are de- 
termined by a matching with the solution in the region of receptivity which is the 
focus of this paper. 

We introduce a Cartesian coordinate system { c*, y*} centered at an arbitrary point 
within this localized region (Fig. 2). Since we are considering the two-dimensional 
case, it seems appropriate to simplify' the formulation by' having the streamfunction 
as the only dependent variable. The streamfunction i j> is normalized by U^S , 
and we define the local coordinates in the streamwise and transverse directions as 
X = x * /6* and Y = y* ft)*, respectively. Although for very large Reynolds numbers, 
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the instability wavelength L* is of the order of several displacement thicknesses, 
the choice of 8* as the local length-scale is motivated by convenience. The slow 
streamwise coordinate representing the growth of the Blasius boundary layer is 
defined as x — x* /('*. The nondimensional time t and frequency are assumed to 
have been normalized by 8* /U^ and its inverse, respectively. 

The streamfunction within the local region, </’ satisfies the two-dimensional 
Navier-Stokes equation, 

<9V 2 0( J) . a0«>avv (i) di^d^_ J_ V7 .i, l(j) _n f01) 

dt + dY dX OY 3 dX R t i. ' ’ ' 

where j=l corresponds to the case of receptivity due to wall suction variation, 
j=2 to the wall admittance problem while j=3 corresponds to the wall hump case. 
Streamfunctions and 0^ satisfy the inhomogeneous wall boundary conditions 


dX 




(2-2a) 


and 

= 'fA ?F'%X) . (2.26) 

corresponding to specified distributions of the mean and unsteady components, 
respectively, of the normal velocity at the wall. In addition, rf.'^ and xj> ^ 2 * must 
satisfy the no-slip boundary condition 


diP U) (X,0) 

dY 


o , i = 1,2 . 


(2.2c, d) 


On the other hand, satisfies the impermeability as well as the no-slip boundary 
condition at the deformed surface location, 


at Y w = ^F*%X) . (2.3a. 6) 

The functions F^\X) (j = 1 — 3) denoting the normalized spatial distributions 
of the wall inhomogeneities are assumed to be arbitrary throughout the analysis. 
Finally, we require that all flow quantities approach tie' respective free-stream values 
far away from the wall. This implies 

V>°') -* (1 + e /s f- iu ") Y as Y -» oo (2.4) 


One may observe that four small parameters appear in the problems defined by 
equations (2.1)-(2.4), namely, R^,ei,e fs and . If one is interested in only a 
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narrow range of the frequency parameter cj, say in vicinity of one of the two neu- 
tral branches, then the two parameters 7?^.' and f/ are related to each other for 
sufficiently large values of the Reynolds number. Similarly, there are restrictions 
on the wall inhomogeneity amplitudes e u v in terms of Reynolds number scaling in 
order that the disturbance produced due to the wall inhomogeneity is a linear per- 
turbation of the base flow. However, in the interest of having a general approach 
for localized receptivity problems, we treat these parameters as independent of each 
other. It will be seen from the following analysis that treating these parameters as 
independent does not affect the amplitude of the generated instability wave. 

Hence, we expand the streamfunction within the local region as 

* (i) «*oC *,K) + e f M-r,Y)c-^ (X t Y) 

+ 0{e}„eW\e h Rj. 1 ) , (2.5) 

where the upper and lower case variables correspond to the steady and unsteady 
terms, respectively. The subscript 0 denotes the base flow quantities, i.e., flow within 
the local region in the absence of any short-scale variation, while the subscript 1 
represents the short-scale perturbations to this base flow due to the local wall in- 
homogeneity. Thus, >P 0 (r,r) corresponds to the Blasius streamfunction, while 
the term involving <!'q(x, Y) represents the forced signature of the unsteady free- 
stream disturbance within the Blasius boundary layer. The quantities T | 1) (.Y, Y) 
and (.Y, K) correspond to the short-scale mean flow perturbations due to local 
variations in the wall-suction and wall-height distributions, respectively. Similarly, 
the term involving represents the short-scale corrections to the unsteady 

signature ^q(x^Y) due to the local inhomogeneity. Specifically, 0^ and 0^ de- 
note the leading order corrections due to interaction of the basic unsteady solution, 

0 0, with the short-scale mean-flow perturbations corresponding to ^ and 
respectively. On the other hand, the motion corresponding to xj?\ ’ ' ls induced di- 
rectly due to the short-scale variation in the unsteady wall-flux. Thus, in each case, 
01 denotes the leading order solution for the short-scale, unsteady motion which 
contains both temporal as well as spatial scales matching those of the instability 
wave. 

The unsteady perturbation 0o to the Blasius boundary layer due to a small ampli- 
tude, time harmonic perturbation to the uniform frec-stream was investigated by 
Lighthffll 11 ), and later by Ackerberg and Phillips^ 12 ! and Goldstein W. Ackerberg and 
Phillips analyzed the unsteady motion in the “leading edge region” corresponding 
to distances of the order of a convective wavelength downstream of the leading edge, 

1. e., .t i = u*(£* + x*)/U£o — 0(1). They showed that for e 1 = (^rV) 1 ^ << 1, 
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the unsteady motion in this region is governed by the linearized unsteady bound- 
ary layer equations. Because of the parabolic nature of these equations, the far 
downstream limit (.iq >> 1) of if.'o corresponds to a particular solution ^Op which 
exhibits a two layer structure in the transverse direction. Acker berg and Phillips 
developed a composite expansion for V’ujj which is valid to 0(.iq) -3 ' 2 in both these 
layers. 



(2.6a) 


where 


«o p = 1 - e 


, 3/2 


2/.I-, 


>}F'M 


3/2 


*Fj},( Ok 

2 5 / 2 j-j /2 


3 3 

( -a -| 7: 

2 2v// 


<J 


+ ;V 3 + oc-5-). 


( 2 . 66 ) 


Here a and >j denote the Stokes layer coordinate s/ojRt,* Y, and the Blasius variable 
(b* s/Rte /C*)Y, respectively. Thus, to the leading order, f/’Op is independent of 
the st teamwise coordinate as well as the mean flow, being identical to the Stokes 
shear wave solution for a purely oscillating flow over an infinite flat plate. Ackerberg 
and Phillips solved the linearized unsteady boundary layer equations numerically to 
demonstrate that the unsteady solution approaches the far downstream behaviour 
through damped oscillations corresponding to the eigensolutions of the linearized 
unsteady boundary layer equations first found by Lam and Rott^ 13 ! . 

Goldstein^ analyzed the unsteady motion in the region farther downstream (aq = 
(){(' 2 )), where the unsteady motion satisfies the Orr-Sommerfeld equation with 
slowly varying coefficients. This is also the region of interest in this paper since we 
assumed the wall inhomogeneity to be located in a region where the quasi-parallel 
stability theory is valid. Goldstein showed that the particular solution (2.6a) is also 
valid in the Orr-Sommerfeld region. However, the Lam and Rott eigensolutions now 
match onto the discrete spectrum of the Orr-Sommerfeld equation, the first mode 
from which corresponds to the Tollmien-Schlichting (T-S) instability wave. This 
is precisely how an instability wave is generated by the leading edge receptivity 
mechanism alluded to in the introduction. Upstream of the location correponding 
to the lower branch neutral stability point, where the instability wave decays expo- 
nentially, the principal contribution to V’o is still provided by the particular solution 
(2.6a). In contrast, if the wall inhomogeneity is located far downstream of the lower 
branch, the contribution from the instability wave generated close to the leading 
edge may become comparable to that from <{>Qp- Since we are interested in the gen- 
eration of the instability wave and not its scattering due to a wall inhomogeneity, 
we ignore any contribution to t/q ( .r , Y ) from eigensolutions generated upstream of 
the wall inhomogeneity. 
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Noting that the base flow quantities \Po and 0 Op depend only upon the global 
coordinate x in the streaimvise direction, we can expand them both as Taylor series 
in x about the origin to obtain 

*o(z,F) = * 0 (F) + #0 AY) x + 0(x 2 ) , (2.7a) 

and 

00 ,(x,Y) = 0op(F) + 0 Opx (F) X - 0(x 2 ) , (2.7 b) 

where subscript x represents the partial derivative with respect to x and the su- 
perscript * denotes the profile of a function at the origin. Since x = O(et) within 
the local region, Eqs. (2.7a,b) imply that the streamwise variations of ^q(x 9 Y) 
and in the local region can be neglected to 0(e/). Thus, both the steady 

and unsteady base flow solutions within the local region correspond to a parallel 
shear flow given by their respective profiles at the origin. This leads to a consider- 
able simplification in the solution procedure for the steady and unsteady short-scale 
perturbations and 0, , j = 1 — 3. The coefficient functions in the governing 
equations for these local perturbations arc now independent of any streamwise coor- 
dinate, thus allowing these equations to be reduced to ordinary differential equations 
after taking a Fourier transform. The rest of this section essentially considers the 
solution for the short-scale perturbations F) and %b[ 2 \X,Y). 

First consider the mean flow perturbations, $^(.Y,F) and v k^(A”, F), due to the 
local variations in wall suction velocity and wall height, respectively. Note that 
the wall admittance variation is manifested only through the unsteady boundary 
condition (2.2b). Therefore, to the leading order, there is no perturbation to the 

/ r> \ 

mean flow for the j = 2 case, and j (A, F) = 0. Substituting the perturbation 
expansion (2.5) into (2.1)-(2.4), collecting terms of O(f^), and transforming the 
resultant equation using the following definition for the Fourier transform, 

OO 

SM =4 = / MX ■ (2.8) 

— OO 

one finds that ^^(F) and 0 j 3) (F) satisfy the time independent form of the Orr- 
Sommerfeld (0-S) equation in the wavenumber space 

~ al)h ~ “ 2) 4,1 = 0 • < 2 - 9 > 


subject to the inhomogeneous set of wall boundary conditions 

^(O) = -F (l) (a)/ia , * ( 1 1) '(0) = 0 , 


(2.10a, b) 
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and 

^ ( , 3) (0) = 0 , ^'(O) - — i'{[(0) F (:t) { o ) , (2.11a, b) 

where the primes denote derivatives with respect to Y and the boundary condition 
for 'I' has been shifted to Y = 0 using a Taylor expansion in Y. Both and 

/ n \ 

vp, ' satisfy homogeneous boundary conditions far aw ly from the wall, 

g ; 0) _ gd 3 )' = o , j — 1,3, as Y — » oo. (2.12a, b) 

Although the term involving j n (2.9) is nominally 0(R ^. 1 ), it has 

been retained in order to satisfy the no-slip boundary condition at the wall. This 
reflects the singular nature of the problem. Equation (2.9) contains additional 
terms which are uniformly 0(RJ» 1 ) and, therefore, cannot be justified on a ratio- 
nal basis. A strictly rational perturbation scheme for the short-scale perturbations 
= 1,3 corresponds to the well-known triple deck theory. However, it involves 
different expansions in three separate regions in the direction normal to the wall 
(Stewartsonf 11 !). In contrast, Eq. (2.9) describes the motion in the entire local re- 
gion and also accounts for the finite Reynolds number effects more accurately than 
just a single term triple deck expansion. 

The time independent Orr-Sommerfeld equation along with boundary conditions 
(2.11a,b) and (2.12a,b) also describes the perturbation to a parallel shear flow over 
an infinite flat plate due to small amplitude, sinusoidal variations in wall height. Ap- 
proximate, analytical solutions to this problem were first obtained by Benjamin^ 15 !, 
He also treated the case of a flow over an isolated bump as a superposition of 
the wavy wall solutions. Subsequently, Lessen and Gangwand 16 ! and Aldoss and 
Reshotkob 7 ! integrated the Orr-Sommerfeld equation numerically to calculate the 
stationary wavy-wall solutions. We also find the solutions for and numer- 
ically by using a spectral collocation scheme on a staggered grid which has been 
described in detail by Macaraeg et atf 18 l. Note that the total mean flow perturbation 
can be obtained by computing the solution for all Fourier wavenumbers and then 
inverting the Fourier transform. However, the focus of this work is to determine the 
amplitude of the generated instability wave. It will be seen that with this limited 
objective, it is sufficient to solve (2.9)-(2. 12) for a single Fourier mode corresponding 
to the (complex) wavenumber of the instability wave in the local region. 

Having obtained the short-scale mean flow perturbation, we now analyze the short- 
scale correction, V’l \ to the unsteady base flow solution. Substituting the pertur- 
bation expansion (2.5) into (2.1), collecting terms ol e/ a ), and introducing 

the Fourier transform (2.8) leads to the 0-S equation with a source term which is 
nonzero for j = 1 and j = 3 

-M ~ ~ “ 2 )d J) + <«*l( 


12 


1 


d 2 


(jb 2 -° 2 ) ^ J) = ~ io Wop(^2 


d 2 


R 6 .d l 




(2.13) 


As seen from (2.13), the forcing term on the right hand side arises due to the inter- 
action of the short-scale mean flow perturbation with the unsteady base flow term 
0 O . Since the wall admittance variation does not cause any mean flow perturbation, 
the source term is identically zero for the j — 2 case. 

Utilizing the perturbation expansion (2.5) along with (2.2)-(2.3), one finds that 0> j 
satisfies the homogeneous boundary conditions at the wall 


0) = V 5 ! ,) '(0) = 0. 


(2.14a, 6) 


Although 0j 2 ^ satisfies a homogeneous governing equation, it lias to satisfy the in- 
homogeneous boundary condition corresponding to the unsteady mass-flux through 
the porous surface 

0< 2) (O) = F (2 ) (o)/i'cv , (2.15a) 

in addition to the no-slip boundary condition 

0j 2), (O) = 0 . (2.156) 

- f Q \ 

On the other hand, shifting the boundary conditions for 0, to Y = 0 leads to 


0J 3) (O) = 0, 0i 3), (O) = — 0o(O) F (3 \k) , 


(2.16a, 6) 


thus implying that 0^ satisfies an inhomogeneous boundary condition as well as 
the inhomogeneous differential equation. The inhomogeneous boundary condition 
(2.1Cb) can be viewed as arising due to the adjustment of the base unsteady solution 

00 to the locally deformed wall geometry. Finally, the streamfunctions 0 j ( j) , j = 

1 — 3 satisfy the homogeneous boundary conditions far away from the wall, 

0^ = 0^' = 0 , j = 1 — 3, as Y — i oc (2.17a, 6) 

Thus, as pointed out by Reshotkot 19 !, the receptivity problem in each case reduces 
to an inhomogeneous boundary value problem vvher ■ the inhomogeneity is char- 
acteristic of the particular receptivity mechanism. C riminale^ 20 ! and TanJ 21 ! also 
developed receptivity theories based on the 0-S equ; tion subject to forcing; how- 
ever, this forcing amounted simply to an inhomogeneous boundary condition cor- 
responding to the free-stream disturbance under consideration and did not involve 
any sliort-scale structure, a necessary ingredient for receptivity as discussed in the 
introduction. On the other hand, the present finite Reynolds number approach is 
similar in spirit to the large Reynolds number asymptotic approach of Goldstein^ 
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and Rubarmh One may, in fact, view the present approach as recasting Goldstein 
and Rulmn’s triple-deck formulation in terms of the more familiar Orr-Sommerfeld 
framework. 

The solution for 0, in the physical space is given by the inverse Fourier integral 

OO 

dV.n=T= j . (2.i8) 

— OO 


where the integration path in the complex k plane is chosen to satisfy the causality 
requirements, i.e., such that the time harmonic solution under consideration here 
corresponds to the time-asymptotic limit of the corresponding initial value problem. 
Since we are only interested in the generated instability wave which propagates 
downstream (i.e., X > 0), the integration contour can be closed in the upper half 
a plane. Applying the Cauchy’s integral theorem shows that the inverse Fourier 
integral corresponds to a sum of contributions from the various singularities of the 
integrand in (2.18). These include the pole contributions from the discrete spectrum 
and branch cut contributions from the continuous spectrum of the Orr-Sommerfeld 
operator. The unsteady motion associated with the generated T-S wave is then 
given by the residue contribution corresponding to the pole at the T-S wavenumber 
aT-s •, the location of which is determined numerically by solving the O-S eigenvalue 
problem. The residue contribution is related to the Fourier transform solution 

by 

,(j) /V V'\ s/2tT l Id ;• _ S' .V 




dj± 


U). 


da 


| o = a T-S 


(2.19) 


We evaluated the right hand side of (2.19) by solving the inhomogeneous O-S prob- 
lem using the spectral collocation scheme of Macaraeg' et a^ 18 ^, and computing the 


derivative in the spectral space, 


0 V 


U ) - 


da 


via a central difference approximation. 


In experiments, it is customary to measure the stream wise velocity fluctu- 

ation associated with the T-S wave. Therefore, using (2.19) and noting that the 
inhomogeneous terms in (2.13)-(2.17) are linear in one can express 

in the following form 




(j)* 

l T-S 


(2.20 a) 


where 

C'i'» = ^}F ln („r-s) A ( /V, ) . (2.204) 

and E q (Y\ to, ) denotes the instability wave eigenfunction for the physical quan- 
tity denoted by q. The eigenfunctions for different, flow variables have been normal- 
ized in such a way that maximum of E U (Y ; o>, /?*. ), the eigenfunction corresponding 
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to the streamwise velocity fluctuation, is equal to unity. The factor F^'(qt-s) is 
the transform of the spatial distribution of the wall inhomogeneity, evaluated at 
the complex wavenumber a'p-s(w, Rd*) °f the instability wave. In contrast, the 
function A* depends upon the external- disturbance frequency, ox, and location of 
the wall-inhomogeneity, Rs-, but is independent of the local geometry. We refer to 
it as the “efficiency function” for the particular comb nation of wall-inhomogeneity 
and free-stream disturbance. 

The product C« relates the amplitude of the geneiated instability wave (at the 
location corresponding to the origin) to the amplituc e of the incident free-stream 
disturbance. Hence, C„ may be x'efered to as the “local coupling coefficient”. 
Since we had remarked earlier about the similarity of the present method with 
the asymptotic approach, it should not come as a surprise that the expi'ession 
(2.20b) for the coupling coefficient is identical in form to the asymptotic results of 
Goldstein^ 2 ! and Choudhari* 5 ’. Equation (2.20b) ecl oes a generic result which is 
valid for all types of linear localized receptivity mechanisms. Essentially, it states 
that the effects of local geometry can be decoupled from the coupling coefficient 
in a very simple manner. This enables one to compare different types of localized 
mechanisms solely on the basis of their efficiency functions A^), irrespective of the 
local geometiy. 

One should note that the values of the coupling coefficient C^\ and hence the 
efficiency function A' J ', are inherently dependent upon the physical quantity chosen 
for measuring the instability wave amplitude. Howcvei , having found the streamwise 
velocity fluctuation pioduced by the generated T-S wave, the residue contributions 
for other flow variables can be determined easily by using the T-S eigenfunction for 
this quantity. For instance, if one measures the pressure fluctuation at the wall, 
Pxl*s(Y — 0), instead of the maximum streamwise velocity fluctuation, one can 


write, 


Pr'ls(-Y,r = o,<) = cy « lar - iV_ "' ) , 

(2.21a) 

where 


&/>-#' F^(a T -s) A«'(u,B..), 

(2.216) 

where the efficiency function Ap is related to A. b ■ 


A« = £?„(!- = 0; A'„» 

(2.21c) 


In the following section we present results illustrating the variation of the efficiency 
functions A^ux, /?*. ) and Ap ;) (u;, ) with the frequency parameter and location 

of the wall inhomogeneity. In computing these results, we have used an approxi- 
mation for the base unsteady solution obtained by using just the leading order 
term in (2.Gb). Computations were also performed usi lg all three terms from (2.6b), 
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and yielded visually indistinguishable results in all cases used for comparison. The 
reasons for truncating the expansion for and its implications on the theoretical 
accuracy of the present approach are discussed in Section 4. 

3. RESULTS 

First consider the variation in the efficiency functions with the frequency of the 
external disturbance for a fixed location of the wall inhomogeneity. As pointed out 
in the previous section, the efficiency function has a different value depending on 
the physical quantity used to measure the amplitude of the generated instability 
wave. To begin with, let us analyze the behaviour of the complex valued 

efficiency functions based upon the maximum streamwise velocity fluctuation. In 
Figs. (3a-3c) and (4a-4c) we have plotted the magnitudes and arguments, respec- 
tively of A n\j = 1-3, as functions of the frequency parameter / = u>*v* jU^ 2 , 
with the local Reynolds number R = y/Rc 7 ^ ns a parameter. For the self-similar 
Blasius boundary layer, R is related to Rfy * , the Reynolds number based on the 
local displacement thickness, via the simple relation R = Rs* /l.7'2. Note that the 
local parameters uj and Rs* were used in Section 2 to formulate the problem since 
the analysis presented herein is valid even for non-similar boundary layer profiles, 
However, for the specific case of a Blasius boundary layer, it is more convenient 
to present the results in terms of the parameters / and R. The symbols in Figs. 
(3a-3c) correspond to results obtained using the present (0-S) approach, while 
the lines represent single term asymptotic predictions based on Goldstein^ and 
Choudharif 5 !. Except in obvious cases or when stated explicitly, the same notation 
will be followed in the remaining figures in this section which involve both symbols 
and lines. The first part of this Section focuses only on the results obtained using 
the 0-S approach, while the latter part presents a detailed discussion regarding the 
comparison of these two approaches. 

Three different locations of the wall inhomogeneity have been considered in Figs. 
(3a-3c). The location nearest to the leading edge corresponds to R = 350, and 
probably represents the lower bound on the range of Reynolds numbers in which 
the present approach can be expected to be valid. The next wall inhomogeneity 
location corresponds to R — 700, while the location farthest from the leading edge is 
assumed to be at R = 1050, i.c., in the range of Reynolds numbers where transition 
is expected to set in on a flat plate without suction. In presence of wall suction, the 
onset of transition may be delayed to much larger values of R. However, since the 
T-S waves leading to transition are usually generated much farther upstream of the 
transition location, receptivity locations corresponding to R > 1050 are not likely 
to he important in practice. For reference, the local growth rate of the generated 
T-S wave at each of these locations has been plotted in Fig. 5 as a function of the 
frequency parameter, /. In addition, the lower and upper branch frequencies at 
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the three stations have been marked on the horizontal axis in Figs. 3a-3c by the 
respective symbols, provided these frequencies fall in the range displayed in these 
figures. 

From Figs. (3a-3o), one can observe two distinct trends in the behaviour of |A U J (/)| 
depending on the type of wall inhomogeneity. For the wall suction (j=l) and wall 
admittance (j— 2) cases, the value of / corresponding to the largest magnitude of 
the efficiency function is less than one fourth of the lower branch neutral frequency, 
fib., at both R = 700 and R = 1050. For larger frequencies, both |a! iJ | and |A„ | 
decrease monotonically until / « f u .b., the frequency parameter corresponding to 
the upper branch of the neutral stability curve. On the other hand, the magnitude 
of A! 3) , the efficiency function in the wall hump case (j=3), increases monotonically 
with / until reaching its maximum somewhere in the vicinity of / = (see also 
Fig. 15b). The rate of increase in |AL 3) | is fairly rapid in the range of subcritical 
frequencies (f < fi.b.)-, but rather slow in the unstable range corresponding to fi.b. < 
/ < fu.b.- One may also observe that the peak value of |A„ | is nearly independent of 
the wall inhomogeneity location, being in the range 0.46-0.48 for all three locations 
of the roughness element. Subsequent to the point of its maximum, |A« (/)| was 
found to decrease rapidly, being smaller than 0.1 for most frequencies larger than 
fu.b.- This behavior suggests that at any given location, the efficiency of receptivity 
due to a wall hump is highest for frequencies which are locally unstable. Finally, 
Figs. 4a-4c show that the arguments of the complex valued efficiency functions 
A = 1—3 change rather slowly in the frequency range of interest. 

In Figs. 6a-6c, the magnitudes of Ap , the efficiency functions based upon the wall 
pressure fluctuation produced by the T-S wave are plotted as functions of / for 
the same locations of the wall inhomogeneity as in Figs. 3a-3c. One may observe 
from Fig. 6a that for small values of the frequency parameter, |Ap 3) | increases at a 
nearly constant rate, in a manner similar to that of |A^| in Fig. 3a. However, the 
value of \Ap l) \ continues to increase until f reaches approximately one-half of the 
lower branch neutral frequency, fi.b.- On the other hand, the function | A „ ^ | , which 
exhibits a boundary layer type behaviour in the range f « fi.b., has its maximum 
near / < fi.b./ 4. Furthermore, Fig. 6a shows that jAp !) | remains virtually constant 
following the initial region of constant slope. This behaviour is markedly different 
from that of |Au^| which decreases monotonically in the range of larger frequencies. 

The behaviour of |Ap 2) |, the magnitude of the efficiency function in the wall admit- 
tance case, is quite similar to that of |Ap'^| corresponding to the wall suction case 
discussed in the previous paragraph. The only qualitative difference between |Ap * | 
and |Ap 2) | is that the latter displays a slight overshoot near / fi.b./ 4, and subse- 

quently, a small dip in the unstable range f Lb . < f < fu.b.- Therefore, the previous 
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discusssion regarding the differences between | Aj/ * | and also holds in the case 

of |Ap 2) | and |A« |. Finally, it is apparent from a comparison of Figs. 3c and 6c 
that the efficiency functions Au and Ap possess different trends even in the case 

/■q \ 

of receptivity due to a wall geometry variation (j—3). The function \A P (/) | does 

/ o \ 

not exhibit the plateau behaviour displayed by |A U (f )\ in the range of unstable 
frequencies, //.t. < / < fu.b.- hi fact, |A^ (/)| has a nearly linear behaviour for all 
frequencies of interest (i.e., / < f u .b.)- Thus, one may conclude that irrespective of 
the type of wall inhomogeneity, the frequency response of a given receptivity mech- 
anism is quite different depending on the specific flow quantity used to quantify the 
receptivity. 

It was shown in Section 2 that the boundary value problem governing the short- 
scale unsteady perturbation involves a different kind of inhomogeneity depending 
on the particular agency inducing the short-scale variation. For instance, the recep- 
tivity due to a wall suction variation is entirely due to internal forcing representing 
the interaction of the Stokes wave with the short-scale? mean flow perturbation. In 
contrast, the variation in the wall admittance leads to receptivity via only a forc- 
ing at the boundary corresponding to the short-scale unsteady mass flux through 
the porous surface. Finally, the boundary value problem for the receptivity due 
to a wall hump contains an inhomogeneous term in both the governing differen- 
tial equation as well as in one of wall boundary conditions. As in the wall suction 
case, the source term in the differential equation arises due to an interaction of the 
Stokes wave with the mean flow perturbation, while the inhomogeneous boundary 
condition represents the direct scattering of the Stokes wave by the wall geometry 
variation. Figures 7a and 7b illustrate the nature of the source term in the dif- 
ferential equation governing the receptivity due to wall suction and wall geometry 
variations, respectively. In both figures, we have plotted the forcing function on 
the right hand side of Eq. (2.13) as a function of the transverse coordinate Y. The 
location of the receptivity region has been assumed o be fixed at R = 1050, and 
results have been plotted for five different values of the frequency parameter /. One 
may observe that for each frequency, the shape of the forcing function in Fig. 7a is 
quite different from that in Fig. 7b. The magnitude of the forcing function in the 
wall suction case is maximum at the wall and decreases rapidly away from it. On 
the other hand, the forcing function in the case of receptivity due to wall rough- 
ness has two lobes. The inner lobe is narrower, occupying less than one third of 
the boundary layer displacement thickness for most frequencies within the unstable 
range. The peak corresponding to the inner lobe also represents the maximum of 
the forcing function across the entire boundary layer. The maximum corresponding 
to the outer lobe is nearly one half of the maximum of the inner lobe. Figures 7a 
and 7b also show that in general, the magnitude of the forcing function is signif- 
icant only within distances less than one half of the boundary layer displacement 
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thic kness away from the wall. Thus, it seems reasonable to conclude that the energy 
transfer from the acoustic wave to the T-S wave is localized to the region close to 
the surface. 

In order to assess the effectiveness of the two types of inhomogeneous terms leading 
to receptivity in the case of wall roughness, we have plotted their contributions sep- 
arately in Fig. 8 along with the sum total. The dashed curve in Fig. 8 corresponds 
to contribution to the efficiency function | A« 3 ^| due to the inhomogeneous term in 
the differential equation, while the dotted curve represents the contribution from 
the inhomogeneous boundary condition. The roughness element is assumed to have 
been fixed at a location corresponding to R = 2800 and the different contributions 
to the efficiency function have been plotted as functions of /. Figure 8 suggests that 
at all frequencies, the contribution due to the inhomogeneous boundary condition is 
much larger than that due to the interaction of the Stokes wave with the mean flow 
gradients. In fact, the latter contribution is nearly out of phase with the former at 
all values of the frequency parameter. Therefore, the magnitude of the efficiency 
function | A„ J> | is approximately equal to the difference of the magnitudes of contri- 
butions from the two inhomogeneous terms. The asymptotic theory (Goldstein^ 2 !) 
predicts that for f » ft.b., the mean flow perturbation becomes very small and 
hence the total efficiency function can be well approximated by the contribution 
due to the inhomogeneous boundary condition alone. However, Fig. 8 shows this 
not to be the case. 

Figures (9a-c) illustrate the variation in |a!/Y./ = 1 — 3, for a fixed frequency dis- 
turbance as the location of the wall inhomogeneity, 7?, is varied. For completeness, 
the arguments of the complex valued functions are plotted as functions of R 
in Figs. lOa-lOc. Results are plotted for four different frequency parameters in the 
range / = 20xl0 -6 to / = 35xl0~ 6 . This particular range was selected based 
upon the corresponding values for the total amplification ratio (the N factor) of a 
fixed frequency T-S wave between the two branches of the neutral stability curve. 
As seen from Fig. (6.4) in Mack^ 22 ^, an amplification ratio of e 2 is attained first by 
a T-S wave with frequency corresponding approximately to / = 27xl0 -6 . The ap- 
proximate N factors for the smallest and largest values of the frequency parameter 
in the range shown in Figs. 9a-9c are equal to 7 and 11, respectively. Hence, accord- 
ing to the well-known c ^ criterion (see Smith and Gamberoni^ 2 '^, Van IiigoiJ 21 ^ , as 
well as Bushncll and Malikl 25 !) the frequencies in this range are likely to be most 
important from the point of view of transition. For reference, the two neutral loca- 
tions for each frequency parameter arc indicated by arrows on the horizontal axes 
in Figs. 9a- 9c. The growth rate variation with respect to the Reynolds number is 
shown in Fig. 11. 

Figure 9a shows that the efficiency of receptivity due to wall suction variation is 
highest in the range of locations which are significantly upstream of the lower branch 
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of the neutral stability curve. The maximum value of |AL^| occurs when the region 
of wall suction variation is located approximately two-thirds of the distance from 
the leading edge to the neutral location. At farther downstream locations, |Ai^| 
decreases, but only to reach a plateau close to the location where the streamwise 
growth rate of the T-S wave is at its maximum. The value of |Au^| remains roughly 
constant between the locations corresponding to the maximum growth rate and the 
upper branch of the neutral stability curve. Shortly after the wave begins to decay 
again, |Au^| drops rapidly to very small values. As remarked before, the behaviour 
of the |aL !) | curve beyond the location corresponding to the upper branch is of 
little practical interest. Figure 9b shows that similar to the wall suction case, the 
local efficiency of receptivity due to a wall admittance variation is highest when the 
variation occurs close to the leading edge. Of course, one must realize that the higher 
levels of efficiency close to the leading edge are offset by the larger magnitudes of the 
decay factors corresponding to the exponential decrease in the T-S wave amplitude 
upstream of the lower branch neutral location. One may further observe that |A^| 
decreases rapidly as the region of wall admittance variation moves away from the 
leading edge. However, the rapid decrease is halted when the wall inhomogeneity 
moves to locations within the unstable region. In fact, the value of | Au^| increases 
slightly for locations close to the upper branch. 

The behaviour of |A^| shown in Fig. 9c as a function of the wall inhomogeneoity 
location is quite similar to that in Fig 3c, where \\ a \ was plotted as a function of 
the frequency parameter with the wall inhomogeneity fixed at a particular location. 
Thus, |Al 3) | increases rapidly until reaching the lower branch neutral location, but 
the increase is quite slow within the range of unstable locations. After crossing the 
upper branch location, |A« | again decreases rapidly to very small values. Note that 
the crude resemblance between the two sets of figures, 3a-3c and 9a-9c, could have 
been anticipated based on the qualitative equivalence between increasing down- 
stream distance for a fixed physical frequency and increasing physical frequency at 
a fixed location. In fact, in the infinite Reynolds number limit, this equivalence can 
be quantified in terms of a similarity parameter R A ! 2 F in the vicinity of branch I 
of the neutral stability curve. 

The present analysis was based on the assumption that the local instability wave- 
length is much shorter than the length scale over which the mean boundary layer 
properties change significantly, i.e., the distance of the wall inhomogeneity from 
the leading edge. This assumption enables one to neglect the variation of the base 
flow over the length of the receptivity region. The validity of this approximation 
can be verified a posteriori from the results in Figs. 9a- 9c. Essentially, the “local” 
approximation is valid if the efficiency function does not vary significantly over an 
instability wavelength. Figures 9a-9c show that the only region within the entire 
domain of interest where the efficiency functions, | A ( ,/ * | , j = 1—3, have a large gra- 
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client corresponds to locations close to the leading edge. For instance, at a location 

* / o ^ 

corresponding to R = 500, the variation in |A„ | over an instability wavelength is 
0(10%) for / = 30e — 6. Consequently, the receptivity predictions in the vicinity 
of this point can be assumed to be accurate only within a comparable bound for 
error. 

In the remaining part of this paper, we compare the results obtained using the 
present approach with predictions based on the asymptotic (triple deck) framework 
used by previous investigators. As remarked before, the coupling coefficient in both 
cases can be expressed as a product of the geometry and efficiency functions, F ( ^ 
and A* J '\ respectively. This allows us to compare the results of the two approaches 
by comparing separately their respective predictions of F (j * and a[/\ For a speci- 
fied spatial destribut ion of the wall inhomogeneity, the geometry factor F^Hax.s.) 
depends only upon the instability wavenumber cxx.S. which is determined by an 
eigenvalue calculation. We already presented results for the imaginary part of ax.s. 
in Fig. 5. This quantity is usually much smaller than the real part of ax.s.- How- 
ever, it may have a significant impact on the value of F 1- J> for certain types of 
geometries. Figure 5 shows that even though the triple deck expansions are valid 
only in the vicinity of the lower branch of the neutral stability curve, the asymptotic 
predictions for the streamwise growth rate, Imag(ax.s.), exhibit correct trends in 
most parts of the unstable region. Of course, since only a leading order solution 
was used to compute the asymptotic results, there are significant quantitative errors 
in comparison to the O-S results. For example, the asymptotic estimates for the 
maximum growth rate are consistently larger than the values predicted using the 
O-S equation. In Fig. 12 we show a comparison of the real parts of the instability 
wavenumber ax—s calculated from the two methods. In contrast to the predictions 
for the imaginary part of cvt-s? one finds that the asymptotic theory always un- 
derestimates the value of Real(ax—S )• The differences in Real(ax-s) niay not be 
crucial for smooth geometries such that the geometry function F (j) (o ) varies slowly 
enough with o. However, for shapes such as a rectangular hump, where F^ J * (a ) 
oscillates on the scale of the hump length, the asymptotic predictions may have a 
significant error. The errors in the asymptotic predictions for both Real(ax-s) an d 
Imag(ax-s) become smaller for larger values of the Reynolds number, especially 
in the vicinity of (Figs. 13a, 13b). However, the extent of this improvement is 
quite marginal even for R as large as 2800. 

Choudliari^ utilized the triple deck theory to examine the influence of the width 
of a suction strip on the. receptivity due to the interaction of a free-stream acoustic 
wave with the mean flow disturbance induced by a uniform suction applied through 
the strip. As discussed above, the local geometry enters only through the factor F ^ 
which is simply the Fourier coefficient of the wall suet ion distribution corresponding 
to the T-S wavenumber, ax-s, Hie specified acoustic frequency, u>. Choudhari^ 
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presented a graph illustrating the variation in 1 ^ with respect to the width of 
the suction strip for four different locations of the strip. The geometry function 
p(D 

was evaluated using the asymptotic value of Q'r-s* However, since the results 
were plotted with the suction-strip width normalized by the local T-S wavelength 
(i.e., as 2nw 3 t T ip/ Rcal((\T-.s))> the same plot can also be utilized in the context 
of the present finite Reynolds number approach by using the value of Real(otT-s) 
obtained from the 0-S eigenvalue problem. The inferences drawn in this manner 
will be exact for the case of a wall inhomogeneity located at one of the neutral 
stations, but only qualitatively true for other locations for which Imcig(aT-s) is 
nonzero. 

In addition to the instability wavenumber o p-S, it is useful to compare the finite 
Reynolds number and asymptotic predictions for E P (Y = 0) which also depends 
only upon the solution to the eigenvalue problem. This quantity represents the 
ratio of the wall pressure fluctuation and the maximum streamwise velocity fluc- 
tuation associated with the T-S wave, and relates the the two efficiency functions 
Aj, j) and i\ ( u }) as shown in (2.23c). Note that in the iisyinptotic theory, one has to 
deal with separate expansions in three different regions in the transverse direction, 
and the maximum of the streamwise velocity fluctuation can only be determined 
by forming a composite expansion which is uniformly valid in all three regions. For 
convenience, the predictions in Goldstein^ and Chouclhari^ utilized an approx- 
imate value for this quantity based upon the maximum across the main part of 
the boundary layer, i.e., across the middle deck of the asymptotic expansion. The 
same approximation has also been used in the asymptotic results used in this pa- 
per. The fact that the maximum of the streamwise velocity eigenfunction across 
the lower deck is slightly larger than this approximate value, especially for fre- 
quencies much larger than the lower branch neutral frequency would indicate that 
the magnitude of E p {Y = 0) based on a composite asymptotic expansion will be 
somewhat smaller than the approximate result displayed here. In Fig. 14 we have 
plotted \E V (Y = 0)| as a function of /, with the wall inliomogeneity location as a 
parameter. One may observe that the agreement between the asymptotic and O-S 
results is quite satisfactory for / < fi.b- However, for frequencies corresponding to 
/ > fl.b.i the asymptotic predictions for E r (Y = 0) are somewhat smaller than the 
values obtained from the O-S equation. Recomputing the asymptotic result using 
the maximum streamwise velocity fluctuation across the lower deck led to reduced 
differences with the O-S predictions in the range / > fi.b.- However, since the cor- 
rect approach involves using the somewhat cumbersome composite expansion across 
all decks, all asymptotic predictions used in this paper are based on the maximum 
across the main deck. 

We now return to the Figs. 3a-3c as well as Figs. 4a-4c, and compare the values 
of the efficiency fuentions A obtained from the asymptotic and finite Reynolds 
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number approaches. Figure 3a shows that for the j = 1 case, i.e., for receptivity 
due to a wall suction variation, the agreement between the two approaches is quite 
good for frequencies close to //.&., the frequency corresponding to the lower branch 
of the neutral stability curve. The agreement is not equally satisfactory in regions 
away from / = fib., but can still be considered reasonable. Figure 15a shows 
that a better overall comparison is possible at larger Reynolds numbers. On the 
other hand, it can be seen from Fig. 3b that in the case of receptivity due to 
a wall admittance variation, i.e., j = 2, the asymptotic predictions for |A„ | are 
in excellent agreement with the finite Reynolds number results even at Reynolds 
numbers as low as 350. The only range of frequencies where the two predictions 
differ significantly corresponds to / << fi.b.- Figure 3c shows that for j = 3, 
viz., receptivity due to a roughness element on the wall, the magnitudes of the 
A ( u 3) function calculated using the two approaches match well at small frequencies, 
/ < fi.b.', however, as seen from Fig. 4c, the arguments of A^ 3 * differ significantly. 
On the other hand, as / becomes larger than //.*.. the two predictions for the 
magnitude of A^ begin to deviate from each other, but the agreement in the phase 
of A ( „ 3) improves substantially. The average difference in the magnitude predictions 
is of the order of 15-20 percent for fi.b. < f < fu.b. Figure 15b shows that this 
difference persists at larger values of tlu; Reynolds number. One may remark at 
this point that the asymptotic predictions of Goldstein^ matched quite well with 
Aizin and Polykov’s^ 6 ! experimental findings related to the receptivity due to the 
interaction of a free-stream acoustic wave with a hump on the wall. However, the 
overall asymptotic predictions were somewhat on the lower side of the range of T-S 
wave amplitudes measured in the experiments. 

The above trends in the comparison of the asymptotic and finite Reynolds number 
approaches would suggest that predictions from the latter approach would match the 
experimental results even better. Finally, if one compares the efficiency functions 
based on the pressure fluctuation at the wall (Fig. 6c), an excellent agreement is 
found at all frequencies as well as Reynolds numbers. The differing comparison for 
the two types of efficiency functions, A ( u j) and Ap j) , can easily be explained on the 
basis of the comparison of E P (Y = 0) discussed previously in the context of Fig. 
14. 

The asymptotic theory has the advantage that it furnishes analytical solutions, 
thus enabling exhaustive parametric studies to be conducted without much cost. 
Therefore, in the absence of detailed numerical results, approximate results can 
be obtained by using the asymptotic predictions. The comparison between the 
two approaches in the preceding paragraphs showed that the triple deck theory 
predicts the overall trends correctly in most of the parameter space. However, 
there are significant quantitative errors in the asymptotic predictions related to 
the neutral curve as well as streamwise growth rates of the T-S waves. This may 
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affect the values of the geometry factor to a considerable extent. On the other 
hand, the asymptotic results related to the efficiency functions were found to be 
relatively more accurate. Therefore, as a first cut measure, the accuracy of the 
asymptotic predictions can be improved by calculating the geometry factor using 
an 0-S result for the instability wavenumber. To investigate if the accuracy of the 
efficiency functions can be further improved by using a little additional information 
from a numerical data base obtained with a finite Reynolds number approach, the 
results from Fig. 3c for R = 700 and R = 1050 are replotted in Fig. 16 after 
normalizing both the ordinate, |A„ |, and the abscissa, /, by their respective values 
at the lower branch of the neutral stability curve. It is obvious that the significant 
discrepancies in Fig. 3c at larger frequencies have now disappeared, yielding an 
excellent agreement between the two approaches at all frequencies as well as wall 
inhomogeneity locations. This suggests a way of profitably utilizing any available 
asymptotic input without sacrificing the accuracy to any great extent. This type of 
correlation, if applicable, will be especially valuable for boundary layer flows more 
complex than the Blasius boundary layer, since the asymptotic theory may be less 
accurate in such cases. 


4. APPLICATION TO RECEPTIVITY IN DISTRIBUTED REGIONS 
OF SHORT-SCALE VARIATIONS 

The discussion in the previous sections was limited to receptivity occuring in local- 
ized regions of wall inhomogeneities. However, the results obtained therein can also 
be extended to receptivity occuring in distributed regions of short scale variations. 
The solution for the localized case can be utilized as part of a slowly varying Green’s 
function for the problem of distributed receptivity in a manner similar to Tarnt 21 ! 
who used this idea earlier in a related context. According to the stability theory 
for weakly non-parallel shear flows, the uniformly valid solution for the streamwise 
velocity fluctuation sufficiently far downstream of a point source at A' = X s is given 
by 


G u (X,Y-J\X a ) = A 0 (x a )A(x-, x a )E u (Y , x ; f)e 


«[e(.Y)-©(.Yj)-^] 


(4.1) 


where the streamwise origin of the coordinate system has been fixed at some arbi- 
trary point sufficiently far downstream of the leading edge. The initial amplitude, 
Ao(:r^), in (4.1) is provided by the solution to the local problem at X=0 (i.e., the 
local coupling coefficient), with the geometry factor JF'^(ar.s.) set equal to unity 
corresponding to a point source excitation. The amplitude and phase functions, 
A(.r; x a ) and O(.Y) are obtained from the weakly non-parallel stability analysis. 
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Summing over the contributions from the entire source region, one would obtain, 


X 


XtT-S 


-/ 


G.(X,YJ\X,)F(X,)dX, 


(4.2) 


where the function F(X a ) now denotes the streamwi.se distribution of the forcing 
function in the equations describing the 0(ff s e^) perturbation. For instance, in 
the wall admittance problem, it corresponds to the specified distribution of the 
unsteady vertical velocity. Thus given the streamwise distribution of the forcing 
function, the integral in Eq. (4.2) can be evaluated to determine the instability 
amplitude at any given station. 

By differentiating with respect to the streamwise variable, Eq. (4.2) can also be 
converted to a wave amplitude equation similar to that obtained by TanrJ 21 !, 


Dut-s 

OX 


. 00 


ut-s + FG U + 0(e(). 


(4.3) 


Equation (4.3) shows that when the receptivity occurs continuously over a large 
number of instability wavelengths, the local change in the amplitude of the T-S 
wave is due to transfer of energy from the mean flow (in other words, the exponen- 
tial growth or decay of instability waves generated upsi ream of the present location) 
as well as due to local input from the external disturbances. As pointed out by Tam, 
the relative contributions from the two types of inputs depend upon the local am- 
plitude of the instability wave. Large values of the instability wave amplitude would 
imply that the contribution due to the external input is much weaker compared to 
the local amplification rate due to transfer of energy from the mean flow via the 
Reynolds stress distribution. Therefore, for locations close to the upper branch of 
the neutral stability curve, the amplitude distribution curve asymptotes to that of 
a pure T-S wave eigensoluion. 

One may also remark in passing that for a spatially periodic forcing such as that 
induced by wall- wavi ness, or suction strips with uniform spacing, the contributions 
from neighbouring locations to the integral in (4.2) would tend to cancel each other. 
Therefore, the Green’s function integral will be dominated by contributions from a 
small range of locations in the vicinity of the point where the integrand is nearly 
stationary. At this location, the instability phase 0 is maximally synchronized 
(or “tuned”) with the phase of the external forcing F thus leading to a minimal 
cancellation in the surrounding region. In other words, the instability wave and the 
forcing function are closest to being in resonance in this region. 
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5. CONCLUDING REMARKS 

An alternative to the triple deck approach was present'd for the prediction of recep- 
tivity due to the interaction of a fiee-stream acoustic wave with localized regions of 
short-scale variations in surface boundary conditions. The receptivity problem re- 
duces to an inhomogeneous Orr-Sommerfold problem in the Fourier transform space, 
and the amplitude of the generated instability wave cai be determined as the residue 
corresponding to the T-S wavenumber. This approach is non-asymptotic, but it has 
the advantage of being relatively more accurate at lower Reynolds numbers. The 
leading order error term is estimated to be 0(/?^* 3 ^ 4 log R $•) or 0(R 6 ^ 5 log Rs*) 
depending on whether the acoustic frequency scales on the frequency corresponding 
to the lower or upper branch of the neutral stability curve. Thus, at finite Reynolds 
numbers, the accuracy of this approach is limited only clue to the “localness” ap- 
proximation and freezing the profile of the undisturbed mean flow over this local 
region. 

The above estimate for the error has been obtained on the assumption that one has 
used the full three term asymptotic expansion in (2.6b) to obtain the unsteady base 
flow solution, ipo p* Replacing if?o P by its leading order behaviour corresponding to 
the Stokes shear wave increases the theoretical estimate for the error to (^(x^ 1 ). 
For sufficiently large Reynolds numbers, this error is of 0{R 6 ^ A ) for / close to 
which is comparable to the error in the single term triple deck expansion. However, 
the error in the finite Reynolds number approach becomes much smaller in the main 
part of the unstable region, since xj -1 is of 0{R b >^) there. Moreover, we compared 
the efficiency functions based on the Stokes wave approximation for 4 } 0p with those 
obtained using the complete three term asymptotic expansion in (2.6a,b) for the 
case of a wall inhomogeneity located at R. = 1050, and tin two results were found 
to be visually indistinguishable on the scale of Figs. 3a, c. 

Since the T-S wave amplitude is determined numerically, using three terms from 
(2.6b) is no more difficult than just the Stokes wave solution. However, a solu- 
tion such as (2.6b) may not always be available, especially in more realistic flows 
corresponding to non-similar base mean flows. In fact, the only available infor- 
mation about the acoustic disturbance may consist of the local fluctuation in the 
free-stream velocity. Since the region of receptivity has been assumed to be a large 
number of convective wavelengths downstream of the leading edge, it seems reason- 
able to compute 0o/> by assuming the airfoil surface to be locally flat and infinite 
in both directions. Then tpop satisfies the O-S equation subject to the specified 
slip velocity in the free-stream. For incompressible flows the Stokes shear wave 
corresponds to the zero wavenumber (i.e., acoustic) solution to the O-S equation. To 
maintain a consistent way of approximating t/>< J;> in all types of boundary layers, we 
have used the Stokes wave approximation even in the case of the self- similar Blasius 
boundary layer. 
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In addition to predicting the generation of the viscous T S waves, the finite Reynolds 
number approach can be utilized to predict the coupling coefficients for the case 
of predominantly inviscid type instabilities in an adverse pressure gradient bound- 
ary layer where the triple deck framework is not applicable. In fact, with minor 
modifications, the present approach can be extended to such complex flows as three- 
dimensional and high-speed boundary layers. These flows are susceptible to various 
types of instability waves, each with its own asymptotic structure. In order to 
compare the effectiveness of a wall inhomogeneity in exciting the different types of 
instability waves, it becomes necessary to have a general approach which will be 
valid for as many of these instabilities as possible. It is for this class of problems that 
the strength of the approach presented in this paper really becomes apparent. Some 
preliminary results on the receptivity in high-speed and three-dimensional bound- 
ary layers using this approach have been presented by the present authors in Ref. 
[10]. Finally, the comparison of the present results with the asymptotic predictions 
suggests that in unsteady problems governed by the triple deck framework, quite 
satisfactory results can be obtained by using only a single term in the asymptotic 
expansion. In addition, such a comparison can also suggest ways of utilizing the 
asymptotic results with better accuracy, especially in more complicated problems 
such as those outlined above. 
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(a) 


A iegion of short-scale variation in wall suction velocity. 



(1>) A region of short-scale variation in wall admittance. 



(c) A region of short-scale variation in wall geometry. 
Fig. 1 Different types of wall inhomogeneitics. 


Receptivity region y* 



Fig. 2 Sketch of the coordinate system. 
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R-1050 



/ 

(a) Receptivity due to wall suction variation 



(b) Receptivity due to wall admittance variation 



/ 

(c) Receptivity due to wall geometry variation 

Fig. 3 Magnitudes of the efficiency functions A u\j = 1 3 from Eq. (2.20b) as 

functions of the non-dimensional frequency param eter / = > with 

the location of the wall inhomogeneity ( R = \f Re t- = Rs- / 1-72) as a 
parameter. The symbols denote the results from the finite Reynolds number 
approach, while the lines correspond to asymptotic predictions. 
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(a) Receptivity due to wall suction variation 
R-1050 



/ 

(b) Receptivity due to wall admittance variation 



/ 

(c) Receptivity due to wall geometry variation 

Fig. 4 Phase of the complex valued efficiency functions A *\j = 1 - 3 from Eq. 

(2.20b) as functions of the non-dimensional frequency p aram eter / = , 

with the location of the wall inhomogeneity ( R = \fRet- — Rs* / 1.72) as a 
parameter. The symbols denote the results from the finite Reynolds number 
approach, while the lines correspond to asymptotic predictions. 



/ 


Fi6 5 s S amew^i,,r Wth ra ‘f’ as a function of / for the 

same waL mhomogeneity locations as in Figs. 3 and 4. 3 
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R-1050 
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(a) Receptivity due to wall suction variation 


FU1050 

R-700 
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(b) Receptivity due to wall admittance variation 


FU1050 

R=700 
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(c) Receptivity due to wall geometry variation 

Fig. 6 Magnitudes of the efficiency functions Ajf\j = 1-3 from Eq. (2.21b) as 
functions of the non-dimensional frequency param eter / = /U^ 2 , with 

the location of the wall inhomogeneity (/? = V^ ( r = / 1-72) as a 

parameter. The symbols denote the results from the finite Reynolds number 
approach, while the lines correspond to asymptotic predictions. 
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(a) Receptivity due to wall suction variation 



(1)) Receptivity due to wall geometry variation 


Fig. 7 The forcing function on the right hand side of (2.13) as a function of 
the transverse coordinate Y = y* ff>* (R = 1050). 
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Fig. 8 Relative contributions of the source term in the differential equation, 

(2.13), and the inhomogeneous boundary condition, (2.16b), to |A^| 
for the case of receptivity due to wall hump located at R = 2800. 
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(a) Receptivity due to wall suction variation 



(b) Receptivity due to wall admittance variation 


F«35e-6 

F-30e-6 



0 1000 2000 3000 

R 

(c) Receptivity due to wall geometry variation 

Fig. 9 Magnitudes of the efficiency functions A = 1—3 from Eq. (2.20b) as 
functions of the wall inhomogeneity location (R = \ARe/« = / 1 .72) with 

/ = ui'uyUZ? as a parameter. 
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(b) Receptivity due to wall admittance variation 



(c) Receptivity due to wall geometry variation 


Fig. 10 Phase of the complex valued efficiency functions A = 1 - 3 from Eq. 

(2.20b) as functions of the wall inhomogeneity location ( R = \ f Rer = 
Rf /1 .72) with / = as a parameter. 



Imag (a t-s) 



Fig. 11 Streamwise growth rate, - Imag (a T -s)> as a function of R for the 
same values of / as in Figs. 9 and 10. 



R=1050 
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Fig. 12 Real part of the instability wavenumber, Heal {a t_ 5), as a function 
of / in the range R = 350 to R = 1050. 
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Fig. 13a Real part of the instability wavenumber, 7?en/ (« r . s ), as a function 
of / in the range R = 1400 to R = 2800. 
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Fig. 13b Streamwise growth rate, — htmy (<v/_.v), as a function of / in the 
range R = 1400 to R = 2800. 
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Fig. 14 Magnitude of the factor E p (Y = 0) from (2.21c), which relates the 

two types of efficiency functions, and Ap -7 ^, plotted as a function 
f ^ or same locations of the wall inhomogeneity as in Figs. 3a-3c. 




(a) Receptivity due to wall suction variation 



/ 

(1>) Receptivity due to wall geometry variation 

Fig. 15 Magnitudes of the efficiency functions A \!\j = 1,3 from Eq. (2.20b) 
as functions of the non-dimensional frequency parameter / = (v*v* fJJ* 
in the range R = 1400 to R = 2800. The symbols denote the results 
from the finite Reynolds number approach, while the lines correspond 
to asymptotic predictions. 
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ft fl.b. 


Fig. 16 Comparison of the numerical and asymptotic results for receptivity 
due to the wall hump : The normalized efficiency factor 

|Au 3) (/)|/|Ai 3) (//.6.)| as a function of the normalized frequency pa- 
rameter f/fi.i. 
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